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Abstract 

I show that the standard picture of the neutron star core containing coexisting neutron 
and proton superfluids, with the proton component forming a type II superconductor threaded 
by flux tubes, is inconsistent with observations of long-period (~ 1 yr) precession in iso- 
lated pulsars. I conclude that either the two superfluids coexist nowhere in the stellar 
core, or the core is a type I superconductor rather than type II. Either possibility would 
have interesting implications for neutron star cooling and theories of spin jumps (glitches). 

Physical Review Letters, in press 



1 



Neutron superfluidity and proton superconductivity in neutron stars could have a number 
of interesting consequences for observed spin behavior and thermal evolution. Interaction 
of superfluid vorticity with the nuclei of the inner crust or superconducting flux tubes in 
the core could lead to the jumps in spin rate, glitches, seen in many neutron stars The 
specific heat of the stellar interior is determined by the state of the matter, while neutrino 
emission processes which cool a young neutron star are strongly suppressed in the presence 
of hadronic superfluids (see, e.g., P). The properties of condensed hadronic systems are 
also of interest in studies of heavy nuclei near the neutron drip line and light halo nuclei 
0]. The properties of hadronic systems in beta equilibrium is therefore a central problem 
in both nuclear astrophysics and nuclear physics. 

Reliable predictions of the pairing states of the neutron star core are not yet possible 
as they require extrapolation of nucleon-nucleon potentials well above nuclear saturation 
density, ps = 2.8 x 10^^ g cm~^. The current picture of the neutron star interior posits 
that the outer core consists of mostly or ^P2-^F2 superfluid neutrons, with about 5% 
of the mass in tvpe II ^5*0 superconducting protons, normal electrons and fewer muons 

Above a density ~ 1.7 ps, the pairing situation is essentially unknown 

The compelling evidence for precession in isolated pulsars j^, 13, 3| provides new probes 
of the state of a neutron star's exotic interior. The periodic timing behavior of PSR B1828- 
11 and correlated changes in beam profile have been interpreted as due to precession with 
a period of ~ 1 yr and an amplitude of ~ 3° jl^ ll^ The measured precession period 
implies a fractional distortion of the star (in addition to its rotational distortion) of e ~ 10^*^. 
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ligj ]. crust stresses [2^ . 



or a 



This deformation could be sustained by magnetic stresses 
combination of the two. 

The picture of the outer core I will consider is as follows. The neutron fluid rotates by es- 
tablishing a triangular array of quantized vortex lines, parallel to the axis of the angular mo- 
mentum of the superfluid and with an areal density of riy = 2mnf2„/7r/i ~ lO'^P(s)^^ cm~^, 
where m„ is the neutron mass, fi„ is the angular velocity of the superfluid and P is the spin 
period. The average vortex spacing is Z„ = n^^^^ — lO^^P(s) cm. If magnetic flux penetrates 
the superconducting core, it is organized into quantized flux tubes, with an areal density 
n$ = i?/$o ~ 10^^i?i2 cm~^, where B12 = 10^^ B, B is the average core field in Gauss and 
$0 = Txfic/ e = 2 X 10^^ G cm~^ is the flux quantum. The average spacing between flux tubes 
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is = n'^^^'^ ~ 3000i?]~2^^^ fm. The magnetic field in the core of a flux tube is approximately 
the lower critical field for the superconducting transition, Hci{— 10^^ G). Unlike the vortex 
array, which is expected to be nearly rectilinear, the flux tube array is likely to have a very 
complicated and twisted structure j2l|. Hence, the vortices are entangled in the far more 
numerous flux tubes. In this Letter I show that this entanglement restricts the precession 
to be of very high frequency and low amplitude, in conflict with observations. 

A flux tube has a core of normal protons; the radius of this region is of order the proton 
coherence length S,p — 30 fm. Outside the flux tube, the magnetic field falls off exponentially 
over a distance equal to the London length: Ap ^ 80 fm. (Type II superconductivity occurs 
when -\/2Ap > ^p). A vortex has a core of normal neutrons, of characteristic radius the 
neutron coherence length, ~ 10 



a vortex magnetizes the vortex 
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"m. Entrainment of protons in the neutron flow about 
|. The length scale over which the vortex's magnetic 
field decays is A„ ^ 10 fm. These length scales are for typical parameters of the core: a 
superconducting transition temperature of 10^ K, an effective proton mass of half the bare 
mass, a proton mass fraction of 0.05 and a total mass density of 3 x 10^^ g cm~^ {2^. The 
protons do not rotate by establishing an array of vortices, but corotate with the crust and 
electron fluid approximately as a rigid body by adjusting the London current [2^ . 

For the angular momentum of the neutron fluid to change, the vortices must move. 
The flux tubes in which they are entangled, however, impede their motion. As a vortex 
segment approaches a flux tube segment, the total magnetic energy increases (decreases) if 
the vortex and flux tube are alig ned (anti-aligned). The energy per intersection is -Ejnt ~ 5 
MeV, with a range ~ Ap 0, bJ ]. The vortices are effectively pinned against the flux array, 
unless they can push the flux tubes through the star, or cut through them by surmounting 
the numerous energy barriers. The pinning force is Fp = Ei^t/Ap ~ 0.1 MeV fm~^ per 
intersection, corresponding to a force per unit length of vortex of fp 

dyn cm~^. Because the core is highly conductive, vortices can push the flux tubes only very 
Slow, 0. Hence, unless t.e vo^t.^es cut th.ou,t the flu. tt.es, tl. core flux tubes and 
vortices move together, and the crust, which is frozen to the strong field emerging from the 
core, approximately follows the motion of the charged fluid. 

The charged fluid of the core responds to changes in its rotation rate approximately 
as a rigid body if magnetic stresses are sufficient to enforce corotation. In the su- 
perconducting core, magnetic stresses propagate as cyclotron-vortex waves of frequency 
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ujcv = {HciB /AuppY^'^k, where k is the excitation wavenumber j25|. Taking k = n/R, where 
R is the stellar radius, pp = 1.5 x 10^^ g cm~^ and B = 5 x 10^^ G, gives the characteristic 
frequency at which magnetic stresses are communicated through the core, cUcd.o — 10 rad 
s~^ The frequency uJcv,o represents an approximate upper limit to the precession frequency 
of the star whole. 

To calculate the precession dynamics, I assume that pressure gradients and magnetic 
stresses force the crust and the charged fluid to move together, and refer to the charged core 
fluid plus crust, whose spin rate we observe, as the "body", even though it contains ^10% 
of the star's mass. By the arguments given above, the core neutron fluid, which accounts for 
^90% of the stellar mass, has its angular momentum fixed to this body, through pinning 
of the vortices to the flux tubes. I assume that the vortices of the inner crust are not 
pinned to nuclei, and nearly follow the body's rotation axis so that they have a negligible 
effect on the precession dynamics [2^. With these idealizations, let us write the the inertia 
tensor of the charged fluid as the sum of a spherical piece, a centrifugal bulge that follows 
the instantaneous angular velocity and an oblate, biaxial deformation bulge aligned with 
the body's principal axis: 



Here /o,c is the moment of inertia of the charged fluid (plus any components tightly coupled 
to it) when non-rotating and spherical, 6 is the unit tensor, is a unit vector along the 
body's angular velocity flc, i^d is a unit vector along the principal axis of inertia, AIq is 
the increase in oblateness about flc due to rotation and AJ^ is the portion of the body's 
deformation that is frozen in the body. Let the neutron fluid's angular momentum vector 
Ln be perfectly tied to the core flux tube array, so that Ln is fixed with respect to the body. 
The total angular momentum is -L = Ln + Lc, where Lc is the total angular momentum of 
the body. The Euler equations in the body frame are /c ■ f2c + x {Lc + Ln) = 0. Define 
principal axes in the body (a;i, X2, x^), where X3 is along the major principal axis (£3 = n^). 
The principal moments of inertia are Ji = Jq + 2AIq/3 — AI^/S = I2 = 13(1 + e)^^ , where 
e = Ald/h > 0. Let the angle between X3 and Ln be a. If Lc, Ln and flc are all aligned, the 
star is in a state of minimum energy for a given angular momentum and does not precess. 
A likely precessional state is one in which Lc and Ln are perturbed slightly about this 
stationary point. To define angles, let X3, Ln, Lc and Oc all lie in a plane at t = 0, with 6 




(1) 
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the angle between xs and L, and 9' the angle between flc and L (see Fig. 1). Linearizing 
Euler's equations in Qci, ^c2 and a, gives the solutions: 

^ci{t) — A COS ilpt — Ooilc: ^c2{t) — Asmflpt, (2) 

where flp = eflc + L^/ h is the body-frame precession frequency, and, 

-e 



A _ Lj^ 



1^1 



The motion of fic is a circle of angular radius |A|/r2c, about an axis that takes an angle 9o 
with respect to X3 in the Xi — x^ plane, completing one revolution in a time 27i/Qp. Removing 
the pinned component (L„ = 0) gives the familiar result of Qp = eQ^, A/Q^ — 0- Restoring 
the pinned neutron fluid, and taking L„ ~ In^c gives flp ~ (/„//i)Oc — lOflc, independent 
of a and 9. When Qp exceeds u;cv,o: as for PSR B1828-11, the precession frequency is likely 
to be closer to cUcvfi — 10 rad s"-*^, still very high. As long as is pinned to the body, the 
star precesses at high frequency for any finite a or 9. 

First consider a state in which Ln and Lc are both aligned at t = 0, but He is not along 
L. In this case, a = 9, giving 

+ «„^«(l-.|). (4) 

The angular velocity vector of the body takes a tiny circle of (angular) radius ~ << ^ 
about axis in Fig. 1, nearly coincident with L. Since X„ is fixed in the body, it too goes 
around 0. For the wobble angle to be ~ 3° with Ln fixed in the body, a << is obviously 
required. Suppose, for example, that a = 0, that is, is aligned with 0:3. Then 

^ -0(l + ^]=-{9 + 9')^-P, 9o = 0. (5) 



Now the angular velocity of the body takes a circle about X3, with radius /3 ~ ^. Relatively 
large-amplitude precession occurs, but still at very high frequency. For precession of large 
amplitude to occur at low frequency, L„ must be able to closely follow the rotation axis of 
the body, so that the 0^ x term in Euler's equations becomes small compared to flc x Lc- 
For this to happen, the neutron vortices must be able cut through the flux tubes. For /? ~ 1°, 
this is likely. The flow of the neutron superfluid past a vortex pinned against a flux tube 
creates a Magnus force per unit length of vortex at location r of — PnK, x ([f2„ — fij x r) , 
where p„ is the mass density of the neutron superfluid and k is a vector in the direction 



of the neutron vorticity with magnitude /i/2m„ |27]. For simphcity, take k = kx^ and 
= ^nXs and Qn = = ^- At t = 0, when Qn, Oc and L all lie in the xi — X3 
plane, the angular velocity of the body is fic = ^(— sin Pxi + cos/Jxs). The instantaneous 
Magnus force per unit length of vortex as a function of position in the star is, for small 
angles, = —XiPnufipx^. If exceeds fp, the pinning force per unit length on a typical 
vortex, the vortices will cut through the flux tubes that are in their way. This condition 
gives |X3| > fp/pnK^f3. For Q = 16 rad s"^ (PSR B1828-11), the inferred (3 of 3° and a 
density = 3 x 10^^ g cm~^; |/^| exceeds fp = 10^^ dyne cm~^ for |x3| > 2 x lO^^i?, that 
is, the Magnus force will force the vortices through the flux tubes almost everywhere in the 
star. This process is highly dissipative. 

As a vortex is forced through a flux tube, quantized vortex waves, kelvons, are excited, 
which propagate along the vortex and eventually dissipate as heat j28|. In the rest frame of 
a straight vortex along the z axis, suppose a straight flux tube in the y — z plane approaches 
at speed v. Since A„ < Ap, the finite (magnetic) radius of the vortex can be ignored. Let 
the vortex and flux tube overlap at t = 0. The vector separation between a point at the 
center of the flux tube which will coincide with the vortex at t = is s(t) = vtx. As a 
simple model of the interaction force, consider fint(s(t)) = Fp(s/Ap)exp[(l — s"^ / Ap)/2]6{z)x, 
where S{z), the Dirac-delta function, gives the distribution of the interaction force along the 
vortex (justified below). 

The relative velocity between vortices and flux tubes is v ~ RQcP in the initial stage that 
flux tubes cut through vortices; taking /? comparable to the observed wobble angle gives 
V = 10^ cm s^^ for PSR B1828-11. As a flux tube passes through a vortex, it excites kelvons 
of characteristic frequency Uq = v/Ap. Kelvons on a free vortex are circularly polarized 
waves. The frequency of a kelvon is related to its wavenumber by ujk = hk'^/2p where n 
is the effective mass ot a kelvon, gwen by , = ™„/.A. The d.— e. parameter A 
is ^ 0.116 — ln(A;^„) for wavenumbers in the range ^ << k << [29], which is easily 
satisfied for the characteristic wavenumber ko = {2jiujQ/fiY^'^ of interest. For v = 10^ cm s~^ 
and C,n = 10 fm, these relationships give ko = 5 x 10^^ fm and /i = 0.06m„. 

The total energy transferred to a vortex per scattering in a potential is given in first-order 
perturbation theory by ||28l | 

AE= ^ / dke 

Anpn^P J -00 



dze 



-ikz 



dte'-^'Uiz,t) 



(6) 
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Here f+{z,t) = fint(2, ■ A* where A* = (x + iy)/\/2 is the right- circularly polarized unit 
vector for a kelvon. Since k^Ap << 1, the flux tube exerts a highly localized force along the 
length of the vortex, justifying the use of the (5-function in model for the force. To estimate 



between intersections of a vortex line with a flux tube. The total number of vortices in the 
core is = 27ri?^f2„/ k. Taking a typical vortex length of R, gives a total dissipation rate 
in the core of Nn^RvAE, or 



a lower limit, since the excitation of flux tubes, which is also dissipative, was ignored. 
Different choices for the dependence of the interaction force on s give the same scaling 
on the parameters appearing in eq. ((Zj), with slightly different numerical factors. The 
use of eq. © assumes that cuttings of the vortex at different locations can be treated as 
separate events, with the excitations due to different cuttings adding incoherently. This 
will be the case as long as /cq^* ^ 1- For ^ = 10® cm s^^ and B12 = 1, kol^ is ~ 2, 
so the approximation of kelvons as distinct wavepackets is a somewhat crude one, but the 
lower limit in eq. (jZj) should be a reasonable estimate for the velocities of interest. Taking 
V = 10® cm s~^, fi = 0.06m„, R = 10 km Ap = 80 fm and Fp = 0.1 MeV fm"^ gives a 
dissipation rate of dE/dt ~ 10^^ erg s~^. Now consider the excess rotational energy of the 
precessing star. The energy in the body E^-ot is related to the energy in the inertial frame 
Eq by -Erot = Eq — L ■ Q,. Most of the angular momentum is in the neutrons, so -L ~ Ln- 
The excess rotational energy is thus AErot ^ In^lf3'^/2 ~ 2 X 10"^^ erg. The characteristic 
damping time is = AEj-otidE/dt)^^'^! hr. Over this short timescale, the precession 
damps to small amplitude. When P is ^0.06°, the Magnus force cannot drive the vortices 
though the flux tubes anywhere in the star; Ln is now flxed in the body, and therefore cannot 
follow the total angular momentum, so the star precesses at frequency flp ~ uicv,o — 10 rad 
s~^. In general, then, long-period precession is not possible. 

To summarize, these estimates show that a neutron star core containing coexisting 
neutron vortices and proton flux tubes cannot precess with a period of ~ 1 yr. Since 
Qp = eQc + Ln/Ii, the fraction of the neutron component's moment of inertia that is pinned 
against flux tubes must be ^ e ^ 10^^. Hence, observations require that neutron vortices 
and proton flux tubes coexist nowhere in the star. Either the star's magnetic fleld does not 



the total dissipation rate in the core, take = n. 



1/2 



{^o/By/"^ as the average distance 




(7) 
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penetrate any part of the core that is a type II superconductor, which seems highly unlikely, 
or at least one of the hadronic fluids is not superfluid. This latter possibility appears unlikely 
in the face of pairing calculations which predict coexisting neutron and proton superfiuids 
in the outer core fl E , 3 fl H III | • 

If the core is a type I superconductor, at least in those regions containing vortices, the 
magnetic flux could exist in macroscopic normal regions that surround superconducting re- 
gions that carry no flux. In this case, the magnetic field would not represent the impediment 
to the motion of vortices that flux tubes do, and the star could precess with a long period. 
Perhaps PSR B1828-11 and other precession candidates are giving us the first clue that neu- 
tron stars contain a type I superconductor. Another, strange possibility, is that "neutron 
stars" are in fact composed of strange quark matter 30 1. 

The possibilities discussed above have interesting implications for models of neutron star 
spin and thermal evolution. Glitch models that rely on vortex-flux tube interactions, e.g., 
21 1, would no longer apply, leaving the inner crust superfluid as a possible origin of glitches 
82| |. The URCA reactions, which are strongly suppressed in regions where both neutrons 
and protons are superfluid, could be significantly increased if macroscopic regions of the core 
are normal, affecting the thermal evolution of young neutron stars. 

I thank R. I. Epstein, C. Thompson and I. Wasserman for valuable discussions. This 
work was supported by the National Science Foundation under Grant No. AST-0098728. 
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FIG. 1: The angles defined in the text, fie takes a circular path about axis o, the dashed line. 
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